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ABSTRACT 

The very existence of more than a dozen of high-redshift (z>4) blazars indicates that a much 
larger population of misaligned powerful jetted AGN was already in place when the Uni¬ 
verse was <1.5 Gyr old. Such parent population proved to be very elusive, and escaped direct 
detection in radio surveys so far. High redshift blazars themselves seem to be failing in pro¬ 
ducing extended radio-lobes, raising questions about the connection between such class and 
the vaster population of radio-galaxies. We show that the interaction of the jet electrons with 
the intense cosmic microwave background (CMB) radiation explains the lack of extended ra¬ 
dio emission in high redshift blazars and in their parent population, helping to explain the 
apparently missing misaligned counterparts of high redshift blazars. On the other hand, the 
emission from the more compact and more magnetised hot spots are less affected by the en¬ 
hanced CMB energy density. By modelling the spectral energy distribution of blazar lobes 
and hot spots we find that most of them should be detectable by low frequency deep radio ob¬ 
servations, e.g., by LOw-Frequency ARray for radio astronomy (LOFAR) and by relatively 
deep X-ray observations with good angular resolution, e.g., by the Chandra satellite. At high 
redshifts, the emission of a misaligned relativistic jet, being de-beamed, is missed by current 
large sky area surveys. The isotropic flux produced in the hot spots can be below ~ 1 mJy 
and the isotropic lobe radio emission is quenched by the CMB cooling. Consequently, even 
sources with very powerful jets can go undetected in current radio surveys, and misclassified 
as radio-quiet AGNs. 

Key words: BL Lacertae objects: general — quasars: general — radiation mechanisms: non- 
thermal — gamma-rays: theory — X-rays: general 


1 INTRODUCTION 

Relativistic jets from powerful radio-loud Active Galactic Nuclei 
(AGN) carry energy and particles from within the sphere of in¬ 
fluence of the central super-massive black hole out to distances 
that far exceed the size of their host galaxy (e.g. Schoenmakers et 
al. 2000; de Vries, Becker & White 2006). When the jet interacts 
with the external medium, a hot spot is formed, powering the ex¬ 
tended structures we call “lobes”. The lobes are characterized by 
relatively low magnetic fields (tens of pG), and are thus inefficient 
synchrotron radiators. Minimum energy arguments, based on the 
assumption of equipartition between the particle and magnetic en- 
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ergies, suggest that the lobes of the most powerful sources can store 
up to 10 61 erg in energy. 

Even when we do not have (yet) resolved radio maps of the hot 
spots and of the lobes (e.g. when we have only single dish obser¬ 
vations), their presence can be inferred from steep low frequency 
radio spectra, characteristic of optically thin synchrotron emission, 
in contrast to the flat spectra of partially self-absorbed, optically 
thick jets 

Fig- 0 illustrates the point by showing the broad band spec¬ 
tral energy distribution (SED) of 1224+2122 (a.k.a. 4C +21.35), 
a 3=0.434 blazar detected in the 7-ray band by the Fermi satellite 
(Ackermann et al. 2011; Shaw et al. 2012)0. Multiple data points at 

1 All shown data points are taken from the ASI Science Data Center 
(ASDC) archive. 
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Figure 1 . SED of CRATES J1224+2122 (alias PKS 1222+21 alias 4C 
+21.35), a FermilLAT detected blazar (Shaw et al. 2012). The solid red line 
is the emission from the single-zone jet model we adopted (see text and 
Ghisellini & Tavecchio 2015), the short-dashed black line is the contribu¬ 
tion from the thermal components (i.e., accretion disc, the torus and the X— 
ray corona), while the solid green line is the sum of the two. The adopted jet 
model assumes that most of the emission comes from a compact jet region, 
and therefore the synchrotron emission is self-absorbed up to ~ 10 3 GHz. 
At lower frequencies the flux is produced by several larger jet regions, self¬ 
absorbing at smaller frequencies, and by the extended radio-lobes, whose 
synchrotron emission is optically thin. The latter component is relatively 
steep (Tj, oc u~ 0 7 ), while the partially opaque emission from the jet is flat 
(F„ oc u°), as illustrated by the long-dashed green lines. 


the same frequency indicate the strong variability typical of blazars, 
that can reach an amplitude of 1-2 orders of magnitude at X-ray 
and 7 -ray energies, while it is much less extreme in the radio band. 
The signature of isotropic and optically thin synchrotron emission 
produced by the hot spots and by the radio lobes is the relatively 
steep slope of the radio spectrum below a few GHz, F v oc zz~ 0 ' 7 - 
zz -1 ’ 0 . Above a few GHz the radio spectrum, produced by the par¬ 
tially opaque jet (labelled as “extended jet” in Fig. |T}, is instead 
much flatter, i.e., F v oc u°. Note that the long-dashed lines at the 
lowest radio frequencies are not fits to the lobe and extended jet 
emission. Their purpose here is simply to guide the eye. Note that 
the case of 1224+2122 is not common, since most blazars do not 
show the spectral steepening at low radio frequencies, even when 
going down to 74 MHz, as shown by Massaro et al. (2013). 

Whereas the core jet emission is strongly collimated and 
boosted by relativistic beaming, the hot spot and lobe emission is 
largely isotropic (or at most mildly beamed). This implies that, as 
a first approximation, the jet-to-lobe flux ratio (i.e., the so-called 
“core dominance”) is an indicator of the observer’s viewing an¬ 
gle. For each source whose jet is pointing at an angle 9 V ^ 1/F 
(where F is the jet bulk Lorentz factor), and whose power is large 
enough to produce radio lobes, there must exist ~ 2 F 2 sources with 
strongly misaligned (and thus undetectable) jets, and whose radio 
lobes should be visible regardless of the viewing angle, given the 
required instrument sensitivity. If so, then radio-galaxies and mis¬ 
aligned radio-loud quasars ought to be 2 r 2 times more numerous 
than blazars. 


Volonteri et al. (2011) showed that this straightforward predic¬ 
tion is confirmed up to z ~ 3, beyond which the number of radio¬ 
galaxies dramatically drops. More specifically, the expected num¬ 
ber density of misaligned radio-loud quasars and radio-galaxies as 
inferred from the Swift/ BAT sample of luminous, massive blazars 
(Ajello et al. 2009), overestimates the number of observed lumi¬ 
nous, radio-loud AGNs (from SDSS-DR7+FIRST) by a factor ~ 3 
in the redshift bin z = 3—4, and by a factor >10 between z = 4—5 
(see Table 1 and Fig. 3 in Volonteri et al. 2011; qualitatively sim¬ 
ilar conclusions are reached by Kratzer & Richards 2015; see also 
Haiman, Quataert & Bower 2004 and Me Greer Helfand & White 
2009). 

Volonteri et al. (2011) put forward some possibile explana¬ 
tions of this discrepancy, namely: (i) heavy optical obscuration of 
high-z radio galaxies; (ii) substantial drop in the average F for 
high-z sources; and (iii) substantial dimming of the radio lobes 
at z>3. The first scenario implies the existence of a large popu¬ 
lation of infrared-luminous, radio-loud quasars with weak optical 
counterparts. The second scenario implies a very large density of 
the emitting relativistic electrons (leading to a very large jet kinetic 
power; see Ghisellini & Tavecchio 2010). 

With respect to the third hypothesis, the idea has long been 
entertained of Cosmic Microwave Background (CMB) photons af¬ 
fecting the behavior of jetted AGNs (e.g. Celotti & Fabian 2004). 

In a recent work, Ghisellini et al. (2014a) have explored 
specifically how the interaction between the CMB radiation and 
electrons within the jet-powered lobes affects the appearance of 
jetted AGNs at different redshifts. The major results can be sum¬ 
marized as follows: owing to its (1 + z) 4 dependence, the CMB 
energy density Ucmb starts to dominate over the magnetic energy 
density Ub within the lobes above z ~3, thereby suppressing the 
synchrotron radio flux at higher z (hereafter referred to as CMB 
quenching). At the same time, high-energy electrons will cool ef¬ 
fectively by Inverse Compton losses scattering off CMB photons. 
Combined, these two effects result in a significant enhancement of 
the diffuse X-ray emission - in the form of X-ray lobes - from 
high-z quasars (Celotti & Fabian 2004; Mocz, Fabian & Blundell 
2011). The hot spot radio emission is less affected by the CMB 
quenching, since they are more compact and more magnetised than 
the lobes. However, at low redshifts, their contribution to the total 
radio flux is less than the contribution due to the radio lobes. 

In this paper, we consider all radio-loud sources at z > 4 con¬ 
firmed to be blazars, and construct their broad band SEDs. These 
SEDs have all a radio flat spectrum (where flat means F(v) oc i/~ a 
with a < 0.4). This spectrum is produced by different regions of 
the jet, from the very compact (subparsec) to the extended regions 
(tens of kpc from the black hole, or more). In all cases we see no 
hints of the presence of a steep thin synchrotron component pro¬ 
duced by the hot spots and by the radio lobes at (rest frame) fre¬ 
quencies <1 GHz, hence providing observational support to the 
CMB quenching scenario for high-z jetted AGNs. By applying 
a state-of-the-art jet emission model, we will fully characterise 
the physical properties of the jets. This will allow us to constrain 
the expected properties of both the hot spots and the lobes, whose 
broad band emission can then be modelled under few reasonable 
assumptions (such as energy equipartition between electrons and 
magnetic field). Finally we will address the observability of such 
low-surface brightness radio-lobes from the known population of 
z > 4 blazars at (rest-frame) frequencies at or below the GHz 
band, specifically with the LOw Frequency ARray (LOFAR, van 
Haarlem et al., 2013). 
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Name 

[1] 

z 

[2] 

SDSS+FIRST 

[3] 

R 

[4] 

T'bol 

[5] 

Lciv 

[6] 

M vir 

[7] 

7-BLR,45 

[8] 

A/fit 

[9] 

Ad,45 

[10] 

Ref 

[11] 

032444.30-291821.0 

4.630 

N 

4468 





5e9 

200 

Yuan et al. 2006 

052506.18-334305.5 

4.413 

N 

1230 





3e9 

148 

Fabian et al. 2001b 

083946.22+511202.8 

4.390 

Y 

285 

47.5 

45.0 

8.9e9 

8.8 

7e9 

178 

Bassett et al. 2004; Sbarrato et al. 2013a 

090630.75 +693030.8 

5.47 

N 

1000 





3e9 

68 

Romani et al. 2004 

102623.61 +254259.5 

5.304 

Y 

5200 





5e9 

75 

Sbarrato et al. 2012 

102838.80 -084438.6 

4.276 

N 

4073 





4e9 

120 

Yuan et al. 2000 

114657.79 +403708.6 

5.005 

Y 

1700 





4e9 

114 

Ghisellini et al. 2014b 

125359.62^105930.5 

4.460 

N 

4700 





2e9 

42 

Yuan et al. 2006 

130940.70 +573309.9 

4.268 

Y 

133 

47.1 

45.1 

2.1e8 

46.0 

8e9 

57 

This paper 

132512.49 +112329.7 

4.412 

Y 

879 

47.4 

45.4 

2.8e9 

46.4 

3e9 

55 

This paper; Sbarrato et al. 2013a 

142048.01 +120545.9 

4.034 

Y 

1904 

47.1 

45.0 

1.9e9 

8.16 

2e9 

54 

Sbarrato et al. 2015 

143023.7 +420436 

4.715 

Y 

5865 





1.5e9 

135 

Fabian et al. 2001a 

151002.92+570243.3 

4.309 

Y 

13000 

47.1 

44.9 

3.2e8 

6.6 

1.5e9 

59 

Yuan et al. 2006 

171521.25 +214531.8 

4.011 

Y 

30000 





6e8 

5.5 

Gobeille et al. 2014; Wu et al. 2013 

213412.01 -041909.9 

4.346 

Y 

24000 





1.5e9 

101 

Sbarrato et al. 2015 

222032.50 +002537.5 

4.205 

Y 

4521 

46.9 

45.1 

1.4e9 

10 . 

2e9 

45 

Sbarrato et al. 2015 


Table 1 . Our blazar sample, containing all known blazars at £ ^ 4. Col. [1]: right ascension and declination; Col. [2]: redshift; Col. [3]: flag for belonging 
to the SDSS+FIRST survey; Col. [4]: radio-loudness; Col. [5]: logarithm of the bolometric luminosity (from S11) Col. [6]: logarithm of the CIV emission 
line luminosity (from S11); Col. [7]: black hole mass (in solar masses, best estimate from S11); Col. [8]: Broad Line Region luminosity (in units of 10 45 erg 
s -1 ), obtained from CIV; Col. [9]: black hole mass (in solar masses) from disk fitting; Col. [10]: accretion disk luminosity (in units of 10 45 erg s —1 ) from 
disk fitting; Col. [11]: References. 


We adopt a cosmology with = 0.3, Qa = 0.7, and Hq = 
70 km s _1 Mpc" 1 . 


2 THE BLAZAR SAMPLE 

Following Sbarrato et al. (2013a), we consider all 31 spectroscopi¬ 
cally confirmed z > 4 quasars in the SDSS+FIRST sample owning 
a radio-loudness R = Fsatiz/F ° > 100, where Fsghz and 

' 2500A 

F ° are the monochromatic fluxes at 5 GHz and at 2500 A, 

2500A 

respectively (see Shen et al. 2011, thereafter Sll). Six bona-fide 
blazars (see Table Q} are then identified from this initial sample of 
31 sources based on Swift X-Ray Telescope (Sbarrato et al. 2012, 
2015; Ghisellini et al. 2014b; 2015) and NuSTAR follow-up ob¬ 
servations (Sbarrato et al. 2013b). Two other quasars out of the 31 
had archival X-ray data, showing intense and hard X-ray spectra: 
by fitting their multiwavelength SEDs, we can now classify SDSS 
J 130940.70+573309.9 and SDSS J132512.49+112329.7 as blazars, 
thanks to their small viewing angles and high Lorentz factors. 

In addition, we include 8 other blazars at z > 4 (see Table 
0 } serendipitously identified in X-rays. Among these 8 sources we 
note that 213412.01-041909.9 is in the SDSS+FIRST sample, but 
it was not included in the Shen et al. sample because there was no 
spectroscopic follow-up. A spectroscopic redshift of z = 4.346 
was determined by Hook et al. (2002) for this source. The same 
occurs for 143023.7+420436 and SDSS J171521.25+214531.8, 
whose redshifts (2=4.715 and 2=4.011, respectively) were first de¬ 
termined by Hook & McMahon (1998). The latter has a further in¬ 
dication of its orientation given by extended radio emission (Gob- 
eille et al. 2014). 

To summarize, we consider a grand total of 16 blazars with 
confirmed 2 > 4; 3 of them are at 2 > 5, and 11 objects belong to 


the SDSS+FIRST sample. In Tab.Qjwe list the main references for 
each object, while the SEDs are shown in Fig.[3]and Fig.lffl 


3 EMISSION MODELS 

In this section we summarize the main properties of the models for 
the jet, the hot spot and lobe emissions we use. Full details can be 
found in Ghisellini & Tavecchio (2009, for the jet) and in Ghisellini 
et al. (2014a, for the lobe), and interested readers are referred there. 

Modelling the jet emission allows us to estimate the total jet 
power, that we assume to be conserved along the jet and to power 
the lobes on a much larger scale. A clear estimate of the jet power 
is therefore instrumental in modelling the lobe emission. Other im¬ 
portant parameters characterising the lobes, such as the physical 
extension, the average magnetic field, and the particle energy dis¬ 
tribution are, basically, unknown. We assume physically motivated 
values for them, and discuss how our choices impact on the final 
results. 


3.1 Jet emission 

For the jet emission we adopt a simple, one-zone, leptonic model 
(Ghisellini & Tavecchio 2009). The model assumes that most of the 
observed radiation is produced in a single spherical region within 
a conical jet of semi-aperture -0 = 0.1 rad. The spherical region 
is initially located at a distance Rdiss from the black hole, it is ho¬ 
mogeneously filled with a tangled magnetic field B, and it moves 
with velocity /3c, corresponding to a bulk Lorentz factor F, at an 
angle 8 V with respect to the line of sight. The resulting Doppler 
factor is 5 = 1/[F(1 — /3cos# v )]. Relativistic electrons are in¬ 
jected throughout the spherical region at a constant rate Q( 7 ), with 

2 Data are taken from the papers listed in Tab. [I] and from the ASDC 
archive at http : //tools . asdc . asi . it/. 
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total comoving power (i.e., as measured in the jet comoving frame) 

r 7m ax 

Pe = Vm e c / QGhd'y, (1) 

" 7min 

where V is the volume of the emitting region and 7 m i n and 7 max 
are the minimum and maximum injection energies. The electron 
energy distribution Q( 7 ) is taken as a smoothly connected double 
power law, with slopes si and S 2 below and above the break energy 
7 b, respectively: 


< 9 ( 7 ) = Q 0 


(7/7b)~ 


1 + (7/7b) 


-S1+S2 


r -3 -in 
I cm s 


( 2 ) 


In the following, we will assume that si ^ 1. S 2 > 2 and 7 m i n = 1. 
For this choice of S 2 , the exact value of 7 max is not critical, since it 
corresponds to the end of the steep tails of both the synchrotron and 
the inverse Compton spectrum. Having specified the injection term, 
the particle density distribution N(-y) [cm -3 ] is found by solving 
the continuity equation at a time equal to the light crossing time 
(that is also the time-scale for doubling the source size), taking 
into account radiative cooling and electron-positron pair produc¬ 
tion. Additionally, we consider the presence of a standard, opti¬ 
cally thick, geometrically thin accretion disk (Shakura & Sunyaev 
1973), assumed to emit a total luminosity La = 10 45 Ld ,45 erg 
s -1 . The assumption of a moving blob may correspond to the “in¬ 
ternal shock” scenario, in which two shells, initially separated by 
some distance A R, move with different F and collide at a distance 
« r 2 AR from the black hole (Rees 1978; Ghisellini 1999; Spada 
et al. 2001). Calculating the particle distribution at the light cross¬ 
ing time, when the particle injection stops, corresponds to the max¬ 
imum produced flux, and allows to neglect adiabatic losses, and 
changes in the size, magnetic field and particle density of the blob. 

We further assume that a broad line region (BLR), responsible 
for reprocessing 10 % of La, lies at a distance Rblr = lO 17 !/^^ 
cm. Further out, at R = 2.5 x 10 18 L ^ 2 5 cm, a molecular torus 
intercepts and re-emits a fraction (~ 20 ^t 0 %) of La in the infrared 
band. 

We consider the following emission processes within the jet: 
(i) synchrotron; (ii) synchrotron self-Compton (SSC); and (iii) in¬ 
verse Compton of the relativistic electrons scattering off photons 
produced by the disc, the BLR, and the dusty torus (collectively 
dubbed EC, for External Compton). Photon-photon interactions 
and pair production are also accounted for, albeit these processes 
turn out to be unimportant for our blazar sample, which is uni¬ 
formly characterised by a steep 7 -ray spectrum. 

The source size if’Rdi SS is constrained by: (i) the minimum 
variability time-scale t val , as the size must be smaller than ~ 
c£ V ar<V(l + 2); (ii) the inverse Compton to synchrotron luminos¬ 
ity ratio, which is a function of Pdiss (see Ghisellini & Tavecchio 
2009); and (iii) the location of the peak frequency of the Compton 
component, which, in turn, dictates the nature of the primary seed 
photons (i.e., IR from the torus vs. UV from the broad lines). 

The resulting source size is too compact to account for the ob¬ 
served radio emission, since synchrotron radiation is self-absorbed 
up to (observed) frequencies ~ 10 3 GHz. At lower frequencies, the 
observed flat radio spectrum (P„ oc v~ a , with a < 0.4, and of¬ 
ten a ~ 0 ) is understood to be produced by the superposition of 
the emission from several, larger, parts of the jet. We do not aim to 
model this part of the spectrum, and in all SEDs shown in Figs. 1 
through 6 the long-dashed lines corresponding to Li, oc v° or to 
Li, oc u^ 1 ^ 3 simply guide the eye. The outer jet, responsible for the 
radio spectrum, can be safely neglected because the corresponding 
luminosity is only a small fraction of the total, and therefore we 
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[2] 
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[6] 
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2 

27 
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le3 

57.58 
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000 

1 

27 
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le3 
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57.23 

0 

4 

27 
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le3 
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57.95 

0 

2 

90 
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le3 
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0 

2 

9 
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le3 

57.11 

56.11 

0 

2 

2.7 

97 

le3 

56.57 

56.57 

0 

2 

27 

95 

le4 

56.55 

56.55 

0 

2 

27 

98 

le2 

56.58 

56.57 

0 


Table 2. Model parameters for the lobe (top part) and host spot (bottom 
part) emission of 1028-0844 (z=4.276). For all models we assumed si = 
—1, S 2 = 2.7 and 7max = 10®. Col. [1]: dimension of the emitting zone 
in kpc; Col. [2]: power injected in each hot spot and lobe in the form of 
relativistic electrons in units of 10 45 erg s~ 1 ; Col. [3]: magnetic field in 
//Gauss; Col. [4]: break Lorentz factors of the injected electron distribution; 
Col. [5] and [6]: logarithm of the magnetic and electron energy contained 
in the hot spot or in the lobe in units of erg; Col. [7]: figure where the model 
is shown. 


can assume that the jet retains the bulk of its power up to the lobe 
site. The magnetic field of the extended part of the jet, up to the kpc 
size, should be large enough to be unaffected by the CMB quench¬ 
ing: assuming a jet at z ~ 4 with B ~1 G at a distance d ~ 10 18 
cm and B oc d -1 , we have that the the jet magnetic field starts to 
be smaller than £?cmb for d > 0.4/(r/ 10 ) kpc, where Bomb is 
the magnetic field in equipartition with the CMB (as seen in the 
comoving frame of the jet) given in Eq. [4] If we further assume that 
the self absorption frequency of the different jet components scales 
as d -1 , we expect that CMB will reduce the synchrotron flux pro¬ 
duced by the extended part of the jets for v< 250 MHz. This crude 
estimate shows that the jet emission remains unaltered at all but the 
very small frequencies. 

The jet carries power in various forms, all of which can be 
conveniently expressed as energy fluxes as 

P; = 7rP 2 r 2 /3c[/'. (3) 

Here, (7/ represents the energy density (in the jet comoving frame) 
in radiation, [/', corresponding to the radiative jet power P r , and in 
magnetic field, [/ B = B' 2 /(8n), corresponding to the jet Poynting 
flux Pb . Finally, we need to consider the kinetic energy density of 
the electrons, U' e = m e c 2 J j\ T (j)ydj, corresponding to the power 
carried by the electrons, P e , and protons, P p (we are assuming to 
have one cold proton - i.e. whose kinetic energy is due to bulk 
motion only - per emitting electron). The total jet power, Pj e t, is 
the sum of all these components. 
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Figure 2. SED of the hot spot A (the brightest) of Cygnus A (blue circles), 
compared to the total radio emission (green squares). The blue bow tie refers 
to the X—ray flux of the hot spot A as observed by Chandra (Wilson, Young 
& Shopbell 2000). The dashed red and solid blue lines are the corresponding 
SSC models. For simplicity, we have fitted the total emission with the lobe 
model. The bottom lines show how the emission of the lobes and the hot 
spots would appear if the source were at z = 4. 

3.2 Hot spot emission 

Powerful jets deposit their kinetic energy into the intergalactic 
medium through a termination shock, the hot spot, which in turn 
injects the power received by the jet into the extended lobe struc¬ 
ture. To find out the main properties of the hot spot emission we 
will examine first the hot spot properties of the closest powerful 
Fanaroff-Ryley II radio-galaxy, namely Cygnus A. Then we will 
interpret the small, but resolved, radio emission of the distant blazar 
SDSS J171521.25+214531.8 (1715+2145 hereafter) as due to the 
hot spot emission, to check that if the properties of nearby and dis¬ 
tant hot spots are similar for similar jet power. 

The radio luminosity of the hot spots of Cygnus A is a fac¬ 
tor 5-10 smaller than the radio luminosity of its lobes (Wilson, 
Young & Shopbell 2000; see Fig. 0. The detection of the hot¬ 
spots in X-ray by the Chandra satellite allowed Wilson, Young 
& Shopbell (2000) to estimate the magnetic field of the hot spot, 
B ~ 1.5 x 1CP 4 G, a value close to equipartition with the emitting 
electrons. The corresponding magnetic energy density is ~2,000 
times larger than the CMB energy density at the redshift z = 0.056 
of Cygnus A. The magnetic field in equipartition with the CMB en¬ 
ergy density, Bomb is 

Bomb = 3.24 x 10“ 6 V (1 + z) 2 

= 8.1 x 10 -5 T jr—^ G (4) 

where the T factor accounts for the possibility to have a relativisti- 
cally moving source, that would see the CMB energy density am¬ 
plified by a factor T 2 . If the magnetic field of the hot spots remains 
B ~ 10 -4 G at any redshift, the corresponding magnetic energy 
density is larger than the radiation energy density of the CMB up to 
z ~ 4.5. Consequently, below this redshift, the effects of the CMB 
on the radio emission of the hot spots are weak. 



Log v [Hz] 

Figure 3. SED of 1715+2145, together with the models for the different 
components of the emission: (i) the thermal component, namely the accre¬ 
tion disk, the torus and the X-ray corona (short dashed black line); (ii) the 
emission from the compact and relativistically beamed jet (the blue solid 
line is the sum of the jet and the thermal components), accounting for most 
of the emission, but not for the radio flux, that must be produced by more 
extended regions (but still relativistic) of the jet, forming the flat radio spec¬ 
trum; (iii) the emission from the hot spots (solid black line) assumed to have 
a radius of R s po t = 1.5 kpc, as measured by VLA (Gobeille et al. 2014); (iv) 
the lobe emission (solid red line), assuming that the power injected through¬ 
out the lobe in the form of relativistic electrons is P e ,lobe — 0-lPjet- Mag¬ 
netic and electron lobe energies are in equipartition. For illustration, the 
dashed red line shows the lobe emission for the same P e i 0 b e and magnetic 
field, but now neglecting the effects of CMB photons. The long-dashed blue 
line is not a fitting model, but simply a line to guide the eye. See Table [2] 
for the adopted parameters. Data points are from the ASDC archive, the 
bow—tie in X-rays is from Wu et al. (2013). The orange circle corresponds 
to 4 mJy at 8.4 GHz, and is the flux observed by the two hot spots. The 
vertical yellow line flags the position of the Ly a line. Above this line the 
flux is severely absorbed by the intervening matter, and we do not consider 
the corresponding data in the fit. (see for illustration the case of 1026+2542 
in Fig. [8}. 


To model the emission of both the hot spot and the lobe we 
have followed the prescriptions of Ghisellini et al. (2014), that will 
be summarised in the next section. To fit the hot spots of Cygnus 
A we used B = 1.1 x 10 -4 G, an injected power of relativis¬ 
tic electrons P e = 7 x 10 45 erg s -1 and a radius of 2 kpc. The 
corresponding magnetic energy is Eb = erg, a factor 3.8 smaller 
than the electron energy E e = 1.6 x 10 57 erg. For the lobe emis¬ 
sion, we have used B = 2 x 10 -5 G, same P e and a lobe radius 
of 58 kpc. Magnetic and electron energies are the about the same, 
Eb ~ E e = 3 x 10 59 erg. 

Fig. El shows how the SED of the lobe and hot spot emission 
of Cygnus A would look like if the source were at z = 4. The 
shape (and the luminosity) of the hot spot would remain unaltered, 
since the total radiative cooling would be dominated by synchrotron 
losses even at this redshift. Its flux would be fainter only because of 
the larger distance. On the contrary, the lobe radio spectrum would 
become much steeper, and thus fainter than the hot spot emission, 
due to the enhanced inverse Compton radiative cooling, now dom- 
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inating over synchrotron losses. The enhanced inverse Compton 
losses produce a correspondingly enhanced X-ray luminosity. 

Among the blazars of our sample, 1715+2145, at z = 4.011, 
was observed by the VLA by Gobeille et al. (2014), who report the 
detection of the source at the level of 387 mJy at 8.4 GHz. Most of 
this total flux (383 mJy) comes from the unresolved core, while the 
remaining 4 mJy are produced by two almost co-aligned sources of 
sizes (i.e. diameters) 0.5±0.1 and 0.4±0.1 arcsec, corresponding 
to 3.5±0.7 and 2.8±0.7 kpc. We interpret them as being hot spots, 
given the similarity of their sizes with the size of the hot spots ob¬ 
served in Cygnus A. For simplicity, we assume that both hot spots 
have a radius of 1.5 kpc, and that they contribute 2 mJy each to 
the observed flux. The power of each jet is 5j e t = 1.4 x 10 4 ' 
erg s -1 , and we assumed that the power in relativistic electrons in¬ 
jected throughout one hot spot is 10 % of Pjet- To reproduce the 
observed flux with these assumptions, the magnetic field in the hot 
spot must be 190 pG. This is not far from equipartition with the 
emitting electrons, since the corresponding magnetic and electron 
energies are Eb = 6 x 10 56 erg and E e = 1.4 x 10 s ' erg. Fig. [3] 
shows the SED of the source together with the assumed models for 
the jet, the hot spots and the lobes (for the latter see below). The hot 
spot luminosity is a small fraction of the power injected through¬ 
out the hot spot in relativistic electrons: in 1715+2145 the hot spot 
emits about 10 45 erg s _1 , to be compared to P e ,HS = 1.4 x 10 46 
erg s' 1 . As a consequence, we will assume that the power in rel¬ 
ativistic electrons injected throughout the hot spot and the lobe is 
the same (P e , hs = P e ,lobe)- Since the magnetic field of the hot 
spots is large, their radio emission is unaffected by the CMB, and 
remains the same when switching it off. 


3.3 Lobe emission 

Theoretical models of jet-powered radio lobes rely on high- 
resolution numerical magneto-hydrodynamic simulations. A cru¬ 
cial role in dictating lobe formation and evolution is played by the 
environment in which they develop (see, e.g., Hardcastle & Krause 
2014 for state-of-the-art simulations and references). While low- 
2 radio galaxies inhabit well-known astrophysical environments, 
such as poor galaxy clusters, the boundary conditions for the high- 
2 jets we are interested in are poorly known at best. As a conse¬ 
quence, we have to rely on numerical work and/or the observed 
blazar SEDs themselves to estimate some of the fiducial values for 
modelling the lobe emission. 

Following Ghisellini et al. (2014a), we assume a spherical 
emitting region of radius Piobe, which is homogeneously filled 
with a magnetic field with a coherence length-scale A = 10 kpc 
(see Carilli & Taylor 2002; Celotti & Fabian 2004). Relativistic 
electrons are injected into the lobe with a total power P e ,iobe and a 
distribution with the same functional form Q(y) °f Eq.[2](but with 
parameters different from the jet ones). We are assuming that the 
jet power is essentially constant all along its length, since the radi¬ 
ated power is estimated to be ~ 10 % of -Pjet (Nemmen et al. 2012; 
Ghisellini et al. 2014c). The steady-state electron energy distribu¬ 
tion is found by solving the continuity equation at particle crossing 

time /cross — (Plobe/c)(l T Plobe/A). 

Fig. [3] shows the SED of 1715+2145 and the lobe model as¬ 
suming Piobe = O.IPjet = 1.4 x 10 46 erg s _1 , a size Piobe = 50 
kpc, and equipartition between relativistic electrons and magnetic 
field, giving 5=20 p Gauss for the lobe. For comparison, we show 
how the lobe emission would look like by neglecting the interaction 
of electrons with the CMB. The CMB enhances the X-ray emission 


by ~ 2 orders of magnitude, while dimming the radio emission at 
frequencies >100 MHz. The lobe emission at lower v is instead 
not CMB-quenched. This happens because low energy electrons, 
responsible for the radio flux at iz<100 MHz, do not have enough 
time to cool, even considering the extra losses due to the interac¬ 
tions with CMB photons. Their number, obtained from the continu¬ 
ity equation, is not controlled by the cooling rate, rather it simply 
equals the injection rate multiplied by the time elapsed from the 
start of the injection (see Ghisellini et al. 2014a for full details). 

3.3.1 Constraints on the physical parameters of the hot spot and 
lobe 

In this part of the section we discuss the constraints we can put on 
the physical parameters of the hot spot and lobe, and how the pre¬ 
dicted spectrum depends upon the assumed values. To this aim we 
consider one particular blazar in our sample, namely 1028-0844. 
In all tested models the magnetic field is estimated by requiring 
equipartition with the electron total energy. Equipartition is actually 
supported by observations of classical double FR II radio-galaxies 
in the radio and X-ray bands, and appears as an educated guess 
(Wilson et al. 2000; Croston et al. 2004, 2005; Belsole et al. 2004). 

Size — The radius of the hot spots of the nearby radio-galaxy 
Cygnus A and of the distant blazars 1715+2145 is close to 2 kpc for 
both sources. In Fig. [4] (right panel) we show the emission assum¬ 
ing 5 hs =1,2 and 4 kpc, always requiring equipartition between 
the magnetic and electron energy densities. One can see that the hot 
spot emission is not sensitive to the changes in size. The assumption 
5 e hs = O.IPjet, for this particular source, leads to an overestima¬ 
tion of the observed radio spectrum, once the hot spot emission is 
summed to the flat radio emission of outer parts of the jet. The syn¬ 
chrotron luminosity dominates over the inverse Compton one, and 
their sum is a small fraction of P e ,HS- 

At low redshift, radio lobes form in the overdense environment 
appropriate for a galaxy cluster. However, at z>4 clusters have to 
form yet, and we can argue that lobes of radio galaxies develop 
within a medium at (or close to) the mean cosmic density. It turns 
out that at a ~ 4 the mean cosmic density is comparable to the 
typical density of the intra-cluster medium ~ 1 Mpc away from 
the center of a typical formed cluster at z< 1. Because of this coin¬ 
cidence, we then expect that the lobe sizes of high-redshift blazars 
are not very different from those of very powerful radio sources 
residing in low-z virialized clusters (although some additional ef¬ 
fect may cause the apparent lack of hot spots and lobes at distances 
larger than 100 kpc for FR II radio-galaxies at 3>1 found by van 
Velzen, Falcke & Kording 2015). 

The left panel of Fig.[4]shows the predicted lobe emission as¬ 
suming different values of Piobe- Different values of Piobe lead to 
different values of the equipartition magnetic field (see the values 
in Tab. O and hence of the self-absorption synchrotron frequency. 
Therefore, although the power injected in the lobes is fixed, the 
synchrotron luminosity gets smaller for larger Piobe as this corre¬ 
sponds to a smaller magnetic field. The X-ray luminosity is instead 
nearly independent upon 5i 0 b e , because in the fast cooling regime 
(i.e., when most electrons with 7 > 7 b cool radiatively) the emit¬ 
ted power balances the (assumed constant) injected power. Fig. [4] 
indicates that, under our assumptions, 5i 0 b e has to be >25 kpc to 
keep model predictions consistent with the observed SED. 

Injected power in relativistic electrons — We assume that the power 
injected in relativistic electrons throughout the hot spot an the lobe 
lobe is a fraction of the total jet power 5j e t. The value of this frac- 
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Figure 4. Predicted lobe and hot spot emission from the blazar 1028-0844. The left panel shows the lobe emission for R\ 0 be — 100 kpc (short dashed black 
line), i^iobe — 50 kpc (solid red), and R\ 0 b e = 25 kpc (dot-dashed blue). For all cases we have assumed that the jet deposits in relativistic electrons the same 
power, i.e. 0.1Pj et = 2.7 x 10 46 erg s 1 for each lobe and hot spot. For the right panel the hot spot radius Rus = 4 kpc (dashed black line), Rus — 2 kpc 
(solid red), and Rhs — 1 kpc (dot-dashed blue). The magnetic field is found imposing equipartition with the total electron energy, both for the lobe and for 
the hot spot. The green solid line refers to the jet model (as in Fig.[5]an d|6j. See Table [2]f or the adopted parameters. 
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Figure 5. Predicted lobe and hot spot emission from the blazar 1028-0844 as a function of the injected power in relativistic electrons. For both panels 
Pe/Pjet = 1/3 (violet dot-dashed line); 1/10 (dashed black line); 1/30 (dot-dashed blue line) and 1/100 (solid red line). For all cases we have assumed that 
the lobe radius is 50 kpc and the hot spot radius is 2 kpc. The magnetic field is found imposing equipartition with the total electron energy. A larger P e implies 
a larger equipartition magnetic field and this in turn implies that the synchrotron luminosity increases more than the inverse Compton one. See Table[2]for the 
adopted parameters. 


tion is basically unknown, but we can set P e ,HS = Pe ,lobe since 
radiative losses are not important in the hot spot. 

If magnetic field, proton energy and electron energy are in 
equipartition and in addition we account for the pV work of the 
lobe, then ~10% appears a fair guess. On the other hand, recent nu¬ 
merical work by Hardcastle & Krause (2014; see their Fig. 5), sug¬ 


gests a value close to 1% (characteristic of an environment where 
hot protons dominate the global energetics). 

The right panel of Fig.[5]shows how the emission from the hot 
spots changes by changing P e HS- For P e , hs = 10 Pjet the corre¬ 
sponding magnetic field is relatively small, and most electrons do 
not cool while crossing the hot spot. As a consequence the particle 
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Figure 6. Predicted lobe and hot spot emission from the blazar 1028-0844 as a function of the break energy of the injected electrons yp (see Eq. 2). For both 
panels: = 10 4 (long dashed violet line); 7 b = 3 X 10 3 (solid red line); 7 b = 3 X 10 2 (dot-dashed blue line). For all cases we have assumed that the lobe 

radius is 50 kpc and the hot spot radius is 2 kpc. The magnetic field is found imposing equipartition with the total electron energy. See Table[2]for the adopted 
parameters. 


distribution does not steepen because of cooling, and the spectrum 
is hard. As P e , hs increases, also the equipartition magnetic field 
increases (see Table together with the cooling rate, inducing a 
steepening in the particle distribution and in the produced spectrum. 

The left panel of Fig. [5] illustrates the effect of changing 
Pe.iobe while keeping Piobe = 50 kpc and the magnetic field in 
equipartition with the total electron energy. This forces the syn¬ 
chrotron luminosity to increase as approximatively oc P e 2 i obe , since 
a larger P e ,tobe implies a higher electron density and a higher 
value of the equipartition magnetic field. By comparing the dif¬ 
ferent models in Fig. [5] we conclude that, unless i?i obe 3> 50 kpc, 
assuming P e ,i 0 be in excess of 10% Pjet leads to an overestimate of 
the radio emission. 

Break energy of the injected electrons — Fig. [6] illustrates the ef¬ 
fects of varying the electron energy at the break of the energy dis¬ 
tribution (Eq. [2}. An increase of 7 (, corresponds to an increase in 
the power injected in high energy electrons (which cool more effi¬ 
ciently). This translates into a higher bolometric luminosity and a 
reduced total energy. This is true both for the hot spot and the lobe. 
Fligher values of 7 b also correspond to a larger Compton domi¬ 
nance, with the inverse Compton spectrum peaking at larger fre¬ 
quencies. 

We can now summarise our fiducial values for the hot spot and 
the lobe parameters. We will assume all sources have hot spots 
of size Phs =2 kpc, and that P e , hs = P e ,iobe = O.lPjet and 
Pe,HS = P e ,iobe = O.OIPjet- We will assume equipartition be¬ 
tween the magnetic field and the total electron energy, while for the 
injected particle distribution (i.e. Eq.|3 we will take power-law in¬ 
dices si = — 1 and S 2 = 2.7, and break and maximum Lorentz 
factor 7 b = 10 3 and 7 max = 10 6 , respectively. As mentioned, the 
steep value of S 2 makes the exact value of 7 max not critical. 


4 RESULTS 

The jet, hot spot and lobe modelling described in the previous sec¬ 
tions are applied to the broad band SEDs of our sample of high -2 
blazars. In Fig.[7]and Fig.[ 8 ]we show the observed SEDs, along with 
the corresponding jet and hot spot+lobe emission models. Models 
of hot spot and lobe emission take into account the CMB when 
solving for the electron radiative cooling, and for the resulting en¬ 
ergy distribution and emitted spectrum. 

The parameters required to describe the SED from mm to hard 
X-ray bands are listed in Tableland Table [4] while the derived 
black hole masses and accretion disc luminosities are reported in 
Table |T] Quite remarkably, the resulting best-fit parameters of the 
jet are very similar to what found for the much larger sample of 
FermifLAT blazars analysed by Ghisellini & Tavecchio (2015) and 
Ghisellini et al. (2014c), though none of our high -2 blazars has 
been actually detected by Fermi/ LAT at 7 -ray energies. 

Typical sizes for the dissipation region are in the range (3 — 
15) x 10 1 ' cm, which lies within the BLR (except in the case of 
0839+5112, for which Rdiss is slightly larger). The jet magnetic 
field is found to be ~ 1 — 2 G, the bulk jet Lorentz factor T ~ 
13 — 15, and the jet power 46.4 < log(Pj et /ergs _1 ) < 47.9. 

Once that data are fit to the composite jet/disk/torus model, 
we can estimate the physical parameters of the hot spot and lobe. 
Assuming P e /Pj e t = 0.1, the equipartition magnetic field is esti¬ 
mated to be between 82 and and 510 fiG for the hot spots and be¬ 
tween 7 and 41 /xG for the lobes, while the total energies range from 
log(i?B/erg) = log(E e /erg) = 56.1 to 58 for the hot spots and 
log(i?B/erg) = log(E e /erg) = 58.5 to 60 for the lobes. For the 
Pe,iobe/Pjet = 10 -2 case we find Bhs between 16.5 and 169 /xG, 
Biobe between 2.3 and 13 /xG, and log(f?B/erg) = log(i?e/erg) 
between 55 and 57 for the hot spots and between 57.5 and 59 for 
the lobes. 

The calculated hot spot+lobe emission, for P e ,iobe/Pjet = 
0 . 1 , sometimes exceeds the observed data, requiring that less power 
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Name 

[1] 

z 

[2] 

-^diss 

[3] 

Rblr 

[4] 

P' 

^e, jet, 45 

[5] 

B 

[6] 

r 

[7] 

0V 

[8] 

7b 

[9] 

7max 

[10] 

Si 

[11] 

S2 

[12] 

log Fr 
[13] 

log Pjet 

[14] 

0324-2918 

4.630 

1200 

1140 

0.02 

2.4 

13 

3 

150 

2e3 

-1 

2.8 

46.0 

46.9 

0525 -3343 

4.413 

630 

1643 

0.06 

2.9 

15 

3 

50 

3e3 

0 

2.4 

46.1 

47.5 

0839 +5112 

4.390 

1470 

1336 

0.04 

1.2 

13 

4 

90 

4e3 

0 

2.6 

45.8 

47.2 

0906 +6930 

5.47 

630 

822 

0.02 

1.8 

13 

3 

100 

3e3 

0 

2.5 

46.0 

46.8 

1026 +2542 

5.304 

504 

916 

0.01 

2.3 

13 

3 

70 

4e3 

0 

2.6 

45.7 

46.7 

1028 -0844 

4.276 

840 

1095 

0.08 

1.4 

15 

3 

100 

4e3 

0 

2.6 

46.2 

47.4 

1146+4037 

5.005 

900 

1006 

7e-3 

1.4 

13 

3 

130 

3e3 

0 

2.6 

45.6 

46.4 

1253-4059 

4.460 

480 

671 

0.02 

1.8 

13 

3 

100 

3e3 

1 

2.5 

45.5 

47.2 

1309 +5733 

4.268 

1200 

756 

3e-3 

0.5 

12 

3 

20 

3e3 

0 

2.7 

44.6 

46.2 

1325 +1123 

4.412 

900 

739 

9e-3 

1.9 

12 

3 

250 

3e3 

0 

2.6 

45.5 

46.3 

1420 +1205 

4.034 

360 

725 

6e-3 

2.6 

13 

3 

120 

3e3 

1 

2.5 

45.5 

46.7 

1430 +4204 

4.715 

540 

1112 

0.1 

1.7 

14 

3 

40 

3e3 

0 

2.6 

46.8 

47.9 

1510+5702 

4.309 

293 

636 

0.04 

2.6 

15 

3 

60 

4e3 

0.5 

2.6 

46.4 

47.5 

1715 +2145 

4.011 

396 

233 

0.013 

0.65 

15 

3 

220 

3e3 

1.5 

3.2 

45.3 

47.2 

2134-0419 

4.346 

432 

972 

7e-3 

2.9 

13 

3 

70 

4e3 

0 

2.6 

45.5 

46.6 

2220 +0025 

4.205 

360 

671 

3e-3 

2.4 

13 

3 

120 

3e3 

1 

2.2 

45.2 

46.4 


Table 3. Adopted parameters for the jet models shown in Fig. [8] Col. [ 1 ]: name; Col. [2]: redshift; Col. [3]: distance of the dissipation region from the black 
hole, in units of 10 15 cm; Col. [4]: size of the BLR, in units of 10 15 cm; Col. [5]: power injected in the jet in relativistic electrons, calculated in the comoving 
frame, in units of 10 45 erg s —4 ; Col. [6]: magnetic field in G; Col. [7]: bulk Lorentz factor; Col. [8]: viewing angle in degrees; Col. [9] and Col. [10]: break 
amd maximum Lorenz factor of the injected electron distribution; Col. [11] and Col. [12]: slopes of the injected electron distribution; Col. [13]: logarithm of 
jet power in the form of radiation, in erg s 1 ; Col. [14]: logarithm of the total kinetic plus magnetic jet power, in erg s —1 . The values of the powers and the 
energetics refer to one jet. 


is deposited into the hot spots and the lobes. For P e ,iobe/Tjet = 
10 -2 we have consistency with the existing radio data for all 
blazars. 

We conclude that the available low-frequency radio data do 
not exclude the presence of extended emission in our blazars, and 
that CMB quenching is instrumental in reducing the radio lobe 
emission below current detection limits in all sky surveys. 

Our modelling relies on a number of assumptions on the fidu¬ 
cial values of key parameters, and the resulting SEDs are degener¬ 
ate with respect to, e.g., the values of the magnetic field strength 
and of the source size. On the other hand, deeper observations by, 
e.g., LOFAR, should be able to detect the lobe radio emission in 
most sources, even in the conservative case P e ,i 0 be = 0.01Pj e t- 
Should the radio (and/or the X-ray) flux of the lobes be detected, 
the source physical parameters would be determined with much 
higher confidence. While in principle a simple low-frequency ex¬ 
cess of the radio flux above the extrapolation of the flat jet spectrum 
marks the presence of a lobe, radio maps (and X-ray detections 
and, possibly, imaging; see §4.1 below) would allow us to measure 
the magnetic field and the injected power, and hence derive the total 
electron and magnetic energetics with minimum assumptions. 


5 DISCUSSION AND CONCLUSIONS 

Previous works (e.g., Volonteri et al. 2011) found that the number 
density of radio-galaxies agrees with the expected parent popula¬ 
tion of blazars at z<3. At earlier epochs, however, there is a clear, 
strong deficit of radio-galaxies. One of the possible solutions for 
this puzzle is the (1 + z ) 4 enhancement of the CMB radiation en¬ 
ergy density, which quenches the isotropic radio emission of the 
lobes while boosts their X-ray luminosity. 

We have then constructed a sample containing all known 
blazars at 2 > 4, modelling their jet emission to ultimately assess 
the broad band emission of the lobes. While current low frequency 


observations of high-z blazars do not show any evidence of isotrop¬ 
ically, optically thin synchrotron emission from extended lobes, our 
results indicate that such observations are simply not deep enough. 
We can predict that, for a range of plausible parameters, deeper 
pointings (by, e.g., LOFAR or by the expanded VLA) should be 
able to detect the lobe emission. At the same time, deep X-ray ob¬ 
servations with high angular resolution (by, i.e., Chandra) have the 
required sensitivity to detect and resolve the lobes in X-rays, as the 
estimated lobe angular size is large enough to avoid flux contami¬ 
nation from the point-like nuclear emission. 

A measure of the lobe radio and X-ray fluxes, as well as of the 
lobe angular size, will allow us to estimate the physical parameters 
of the source with no ambiguity and without invoking equipartition 
between the particles and the magnetic field. Besides testing the 
validity of the CMB quenching scenario, this will give information 
about the energy in the proton component of the lobes, and possibly 
the coherence length of the magnetic field. 

According to our scenario, if the lobes were indeed large 
(namely, i?i 0 b e > 25 kpc) and in equipartition, then it would be 
difficult to detect them in any radio-galaxy in the radio band at 
high redshifts. The radiogalaxies at z > 4 with detected extended 
radio emission apparently challenge our proposed CMB quenching 
scenario (see e.g. TN0924-2201 at z = 5.19, van Breugel et al. 
1999; 6 C 0140+326 at a = 4.41, Rawlings et al. 1996; 4C63.20 at 
2 = 4.261, Lacy et al. 1994a; PKS 1338-1942 at z = 4.11, De 
Breuck et al. 1999). We defer a detailed study of these sources to a 
paper in preparation, yet we can anticipate that the reason of their 
strong radio emission most probably lies in the fact that what we 
detect are the hot spots, not the lobes, (see e.g. De Breuck et al. 
2010 for an estimate of the lobe sizes). The strong magnetic field 
of their hot spots (of order of hundreds of /rG), makes synchrotron 
losses competitive with the inverse Compton scattering off CMB 
photons. 

Our results have a profound impact on the estimated fraction 
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Figure 7. SED of the blazar 0324-291 . The blue solid line is the model for 
the sum of jet emission and thermal components (see Fig. |TJ. For the hot 
spots and the lobes we show two possible models: one with 0.1Pj e t (thin 
green lines) and the other with 0.01 Pj et (thin red lines). We show separately 
the contribution of the hot spots and the lobes, as labelled, while the thick 
green and red lines are the corresponding sums. The radius of the hot spots 
and the lobes are assumed to be Phs =2 kpc and (ij 0 |, e =50 kpc. The 
straight dashed blue line interpolating the radio spectrum is not a fit, but 
only a guide to the eye. The hatched grey area in the upper right corner 
shows the sensitivity of Fermil LAT (5<r) after 5 years of observing time. 
We also report (orange line) the sensitivity of LOFAR. 


of radio loud high -3 AGNs. The jet emission in blazars is in fact 
easily visible even at high redshifts and is unaffected by the CMB 
quenching since: (i) jet emission is strongly enhanced by beaming, 
and (ii) the radio flux originates in compact regions where magnetic 
energy density dominates over the CMB. The beaming pattern of 
jet emission decreases fast for increasing viewing angles: for T = 
15, the flux at 9 V ~ 15° is dimmed by orders of magnitude 
compared to the one at 9 V = 3°. This implies that misaligned, 
high-z jets are invisible for the sensitivity of current instruments. 

Fig.[9]illustrates the point. We compare the SED of 1026+2542 
(z=5.3) computed in §4 to the corresponding model as observed at 
viewing angle 9 V = 15°. By modelling the hot spot+lobe emis¬ 
sion as detailed in H3.31 (green line in Fig. [9}, we find that the re¬ 
sulting isotropic radio flux would reach the mJy level at frequen¬ 
cies ~300 MFlz, and would be much weaker at higher frequen¬ 
cies. As a result, a source like 1026+2542 whose jet is observed at 
9 V > 15° would be easily detected in current optical survey, thanks 
to its quasi-isotropic accretion disc emission, but would fail to enter 
radio-catalogs like FIRST, with its sensitivity limit of 1 mJy at 1.4 
GHz. Therefore, slightly misaligned high-z jetted sources would 
be classified as radio-quiet. The very term “radio-loud” becomes 
misleading. The idea that some high-z radio-quiet AGNs could 
have extended X-ray emission, revealing the presence of the jet, 
was first put forward by Schwartz (2002) analysing the X-ray data 
of three z ~6 radio-quiet sources observed by Chandra. The short 
exposure of these observations seemed to suggest that the X-ray 
emission was extended. However, the extended emission of one of 
them was then associated to a galaxy (Ivanov 2004), and this was 
confirmed by longer Chandra observations by Schwartz & Virani 
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Figure 8. SED of the z > 4 blazars. As in Fig. [7] but showing only the 
sum of the hot spot and lobe components. For 0839+5112 the black solid 
line shows the SDSS spectrum dereddened for Galactic absorption. 

(2004). Brandt et al. (2002) did not confirm the extended nature of 
the X-ray flux observed in the other two sources. 

Fig. [9] also shows that most of the lobe radiative losses are 
instead in the X-ray band. Electrons emitting above ~1 keV are 
cooling in less than one source crossing time, and therefore are re¬ 
sponsible for releasing almost the entire injected power Pl° he ~ 
O.IPjet = 5 x 10 45 erg s _1 . Lobes are then radio-quiet and X-ray 
loud. Interestingly, in the 2-10 keV X-ray band the extended lobe 
and the nuclear accretion disk corona emissions can be compara¬ 
ble, as shown in Fig. [9] As a linear proper dimension of 100 kpc 
corresponds at z = 5 to an angular size of 16 arcsec, sub-arcmin 
X-ray imaging is required to disentangle the two components. 
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Figure 8. continue. SED of the z > 4 blazars. As in Fig. 0 but showing 
only the sum of the hot spot and lobe components. 
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Figure 8. continue. SED of the z > 4 blazars. As in Fig. [7] but showing 
only the sum of the hot spot and lobe components. 
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Figure 8. continue. SED of the z > 4 blazars. As in Fig. [7] but showing 
only the sum of the hot spot and lobe components. 
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Table 4. Adopted parameters for the hot spot and lobe models shown in Fig. 
[7] and Fig-Ill For each source, the first raw corresponds to P e /P] e t — 0.1, 
while the second raw corresponds to P e /Pj e t = 10“ 2 . Col. [1]: name; 
Col. [2]: logarithm of the power injected throughout the hot spots and lobes 
in relativistic electrons in erg s 1 ; Col. [3]: magnetic field of the hot spot in 
fi G; Col. [4]: logarithm of the energy in magnetic field contained in the hot 
spot, in erg; Col. [5]: magnetic field of the lobe in //G; Col. [6]: logarithm of 
the energy in magnetic field contained in the lobe, in erg; All sources have 
Rhs = 2 kpc, Ri 0 be = 50 kpc, 7 b = 10 3 , 7max = 10 6 , and slopes of 
the injected electron distribution si = —1 and S 2 = 2.7. The values of the 
powers and the energetics refer to one jet and one hot spot and lobe, while 
the lobe flux shown in the figures corresponds to two hot spots and lobes. 
*: for this source the hot spot is detected, so we have used the observed 
fins =1-5 kpc. The magnetic field is found without invoking equipartition, 
and the electron energy of the hot spot is log(E e ,Hs/ er g) = 57.1. 
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Figure 9. Top panel: the SED of the blazar 1026+2542 at z = 5.3 is 
fitted with a model (parameters in Tab. [3} with Y — 15 and 0 V = 3°. as 
in Fig. [8] Here we also show the contribution from the jet (long-dashed 
violet line), the thermal components (disc, torus and corona, as a short- 
dashed black line) and their sum (solid blue line). The solid green and the 
dot-dashed red lines are the estimated lobe+hot spot emission as in Fig. [5] 
(the two models differ from the assumed injected electron energy, equal to 
10% of 1% of -Pjet, respectively). Bottom panel: if the very same source 
is observed at 9 V = 15° (i.e., only slightly off the beaming angle), the 
jet emission (long-dashed red line) is so weak that it would be missed by 
all current surveys, while the isotropic components (disc, torus, corona and 
the lobes) would remain the same. Current optical surveys, as the SDSS, 
would easily detect the source, and classify it as a radio-quiet quasar. The 
lobe + hot spot emission, slightly below the sensitivity of FIRST (1 mJy at 
1.4 GHz), could be detectable by LOFAR (see Fig.[8j; furthermore, pointed 
Chandra X-ray observations could detect and resolve the X-ray emission 
of the lobes (note that 50 kpc at 2 = 5.3 correspond to an angular size of 
8.2 arcsec in the adopted cosmology). These observations could establish 
the jetted nature of the source. 
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